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Abstract. The main instabilities induced by rotation in stellar interiors are 
O^l ' described. We derive from first principles the general equation describing the 

| transport of the angular momentum. The case of the transport of the chemical 

, species is also discussed. As long as the mass loss rates are not too important, 

■ meridional currents, by advecting angular momentum from the inner regions to 

the outer layers, accelerate the stellar surface during the Main Sequence phase. 
A 9 M Q stellar model at solar metallicity with an equatorial velocity at the 
beginning of the core H-burning phase equal to 340 km s" 1 reaches the break- 
up limit during the MS phase. The model with an initial velocity of 290 km s _1 
^ . approaches this limit without reaching it. The models with 290 km s _1 and 340 

' km s _1 predict enhancements of the N/C ratio at the end of the MS phase equal 

' to 2.8 and 3.2 times the initial value respectively. 
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R ecent discussions of the va r ious instabilities i nduce d by rotation may be found 
in lMaeder fc Mevnetl <|200Cf l. iHeger fc Langerl (1200(1 ). and lTalonl (|2004fi . Among 



the most important instabilities are the secular shear instability and the merid 
ional circulation. These instabilities drive the transport of the chemical species 
and of the angular momentum. Let us briefly recall the physical principles un- 
derlying these two instabilities. 

I n a rotating star local radiative equilibrium cannot be achieved (|Von Zeioell 
K<h ; 119241 : feddington 1925 : Ivbg<lll92fift ■ As a result some parts of the star are heated 

while others are cooled. The buoyancy forces then drive a large scale motion, 
called the meridional circulation. Meridional circulation, contraction/expansion 
of the stellar layers, and convection create gradients of the angular velocity inside 
the star. These gradients produce instabilities known as shear instabilities. The 
physical reason for this instability lies in the fact that the minimum energy 
state of a differentially rotating fluid is solid body rotation. The star will tend 
to approach this state by homogenizing the angular velocity by turbulent mixing. 

In a radiative zone, the vertical stable density stratification counteracts 
both the shear and meridional instability. In that respect the //-gradients play 
a key role as a stabilizing agent. These gradients may even, depending on the 
phys ics involved in the model, completely inhibit the mixing ( Mevnet Sz Maederl 



119971 ). There are different methods in the literature for accounting for the effects 
of the //-gradients on the mixing. Some authors choose a parametric approach 
consisting in multiplying the //-gradient by a free parameter, /„, smaller than 
one in order to weaken the stabilizing effect of the //-gradient. The value of 
is chosen in order to enable the stellar models to reproduce the observed surface 
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enrichments ( Heger &: Langerll2000l ) . The very small values chosen (of the order 
of 0.05) well illustrat es the strong i nhibiting effects of th e //-gradients. Other 
methods devised by iMaederl (jl997f ) , iTalon Zahnl (|l997f ) take into account of 
the effects of the strong horizontal turbulence (see below) and of the heat ex- 
changes which can occur during the mixing process. This reduces the vertical 
drawback force and makes easier the development of the instability. In that case, 
the energy available in the shear can always be used for performing some mixing. 
Expressions for the shear diffusion coefficients can then be deduced without the 
need of artificially reducing the inhibiting effect of the //-gradients. 

In the horizontal direction, i.e. on a isobar, in contrast to what happens in 
the vertical direction, the shear instability is not inhibited by a stable density 
stratification and the turbulence can develop without difficulty. Thus as long 
as horizontal gradients are continuously built up by e .g. meridion al circulation, 
a strong horizontal turbulence develops in the star ((ZahnlH9 92). This strong 
horizontal turbulence will also erode very efficiently any horizontal gradient of 
£1. Thus the star can always be considered to be in a state of shellular rotation, 
characterized by constant v alues o f o n isobars. Starting from this a priori, 
but reasonable hypothesis, IZahnl ( 19921 ) proposed a consistent theory of the 
interaction between the shear instability and the meridional circulation. The 
theory can be kept one dimensional thanks to the hypothesis of strong horizontal 
turbulence. 



2. Equations for the transport of the angular momentum 

Let us express the rate of change of the angular momentum, 4j , of the element 
of mass in the volume ABCD represented in Fig. ^ 

dC A/r 

where M is the momentum of the forces acting on the volume element. We 
assume that angular momentum is transported only through advection (by a 
velocity field U) and through turbulent diffusion, which may be different in the 
radial (vertical) and tangential (horizontal) direction. The component of the 
angular momentum aligned with the rotational axis is equal to 1 

pr 2 sin 9d9dfdr r sin 6fl rsin#, 

Mass of ABCD velocity 

distance to axis 

where O = (p. Since the mass of the volume element ABCD does not change, 
the rate of change of the angular momentum can be written 

pr 2 sm8d6dipdr^-(r 2 sm 2 8n) Mr . (1) 



The components perpendicular to the rotational axis cancel each other when the integration 
is performed over tp. 



Rotation and Mixing in Massive Stars 



3 



Due to shear, forces apply on the surfaces of the volume element. The force on 
the surface AB is equal to 

an 2 

r] v rsin#— r sin 6*d6>cko . 

s or v ' 

vertical viscosity v surface AB 

vertical shear 

The component of the momentum of this force along the rotational axis is 

r] v r sin 9-^-d9dip r sin 6^ 

s v ' distance to axis 

force on AB 

The component along the rotational axis of the resultant momentum of the 
forces acting on AB and CD is equal to 

8 (90 

— (7] v r 4 sin 3 9d8dcp— )dr. (2) 
or or 

The force on the surface AC due to the tangential shear is equal to 

dfi 

7]^ rsmO—^Q r sin9d<pdr, 

v ' surface AC 
tangential shear 

where rj^ is the horizontal viscosity. The component along the rotational axis of 
the resultant momentum of the forces acting on AC and BD is equal to 

d dtl 
-— (rj h r 2 sin 3 9drd<p—)rd6. (3) 
rot) o9 

Using Eqs. E [2]and|SJ simplifying by drdOdip, one obtains the equation for the 
transport of the angular momentum 

d d dft d 80, 

pr 2 sm6—(r 2 sm 2 6n) Mr = Q^(v V r 4 sin 3 6— ) + — (r/ h r 2 sin 3 9— ). (4) 

Setting r\ v = pD v and r/h = pDh and dividing the left and right member by 
r 2 sin#, one obtains 

p-(r 2 S m 2 9n) Mr = __(,IV*_) + — elfe {pD h ^9-). (5) 
Now, the left-handside term can be written 

/O^O 2 sin 2 9Vt) Mr = ^(z 97 * 2 sin2 Q^l)M r ~ r 2 sin 2 9Q^\ Mr - 
Using the relation between the Lagrangian and Eulerian derivatives, one has 

p^-(r 2 sin 2 90) Mr = ^r{pr 2 sin 2 9£l) r + U • V(pr 2 sin 2 0O) - r 2 sin 2 0jA Mr . (6) 
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Figure 1. The momentum of the viscosity forces acting on the element 
ABCD is derived in the text and the general form of the equation describing 
the change with time of the angular momentum of this element is deduced. 
The star rotates around the vertical axis with the angular velocity f2; r and 
9 are the radial and colatitude coordinates of point A. 



Using 

^[M,=|| r + U.Vp. 

and the continuity equation 

g?| r = -div(pU), 
one obtains dp/dt|jw r + pdivU = 0, which incorporated in Eq. Ogives 

p^-(r 2 sin 2 9Q) Mr = ^r(pr 2 sin 2 6Q) r + VCUpr 2 sin 2 0fi). 
dt ot 

Developing the divergence in spherical coordinates and using Eq. [51 one fi- 
nally obtains the equation describing the tr ansport of the angular momentum 
( Maeder fc Zahn1ll998l : iMathis fc Zahnll2004l ) 

d 1 d 1 d 

— (pr 2 sin 2 en) r + ^^(pr 4 sin 2 6w r n) + ^^^(pr 2 sin 3 9wgQ) = 

ot r z or r sin ti ot) 

sm 2 e a , _ 4 sa id frk , 3o on. 

where w r = U r + r is the sum of the radial component of the meridional circula- 
tion velocity and the velocity of expansion/contraction, and wg = Ug, where Ug 
is the horizontal c omponent of the meridional circulation velocity. Assuming, 
as in IZahnl ( 19921 1 , that the rotation depends little on latitude due to strong 
horizontal diffusion, we write 



fi(r,0) = O(r) +n(r,0), 
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with Ci <C Cl. The horizontal average Cl is defined as being the angular velocity 
of a shell rotating like a solid body and having the same angular momentum as 
the considered actual shell. Thus 

- _ jnsm 3 9d6 
~ Jsm 3 6d9 ' 

Any vector field whose Laplacian is nul can be decomposed in sp herical harmon- 
ics. Thus, the meridional circulation velocity can be written ( Mathis &: Zahnl 
2004) 

U = J2 0i(r)fl(co8 0)e r + J2 ^( r ) dP ^ COS ^ e , 

l>0 l>0 

where e r and eg are unit vectors along the radial and colatitude directions 
respectively. Mul tiplying Eq. El by sin fldfl and integrating it over 9 from 
to 7T, one obtains (|Maeder Zahnlll998t l 

^y~V)r = g^O^PW - 5r]) + ^ [pD^-} . (8) 

It is interesting to note that only the I = 2 component of the ci rculation is able 
to ady ect a net amount of angular momentum. As explained in ISpieeel Sz Zahnl 
( 1992) the higher order components do not contribute to the vertical transport of 
angular momentum. Note also that the change in radius r of the given mass shell 
is included in Eq.|Hl which is the Eulerian formulation of the angular momentum 
transport equation. In its Lagrangian formulation, the variable r is linked to 
M r through dM r = 4Trr 2 pdr, and the equation for the transport of the angular 
momentum can be written 

p|(r 2 0) Mr = ^- 2 Upr^U 2 {r)) + \£ ( P D v r^) . (9) 
ot 5r z or r z or V or I 

The charact eristic time associated to the transport of Q by the circulation is 
( Zahnlll992T ) 



tn~t KH l—-\ , (10) 

where g s is the gravity at the surface and txH the Kelvin-Helmholtz timescale, 
which is the characteristic timescale for the change of r in hydrostatic models. 
From Eq. 1101 one sees that tQ is a few times txH, which itself is much shorter 
that the Main Sequence lifetime. 

In order to resolve Eq. |5J one needs expressions for U2(r), D v , and D^. In 
the following we indicate the general lines of the reasoning and give the refer- 
ences where more complete mathematical derivations can be found. The veloc- 
ity of meridio nal circulation is derived from the equation of energy conservation 
i Mestelllflfiah 



pT 



— + (re r + U)-V5 



div(xVT) + pe - divF h (11) 
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where S is the entropy per unit mass, x the thermal conductivity, e the rate of 
nuclear energy per unit mass and the flux of thermal energy due to horizontal 
turbulence. All the quantities are expanded linearly around their average on a 
level surface or isobar, using Legendre function of order 2 i-2(cos 9). For instance 

T(P,6) = f(P)+TP 2 (cos 9). 



Then Eq. ^2 is linearized and an expression for U ?. can be deduced 
Using the same method iMaeder fc Zahnl (jl998f ) revised the expression for U 2 
to account for expansion and contraction in non-stationary models. They also 
studied the effects of the /^-gradients (mean molecular weight gradients), of the 
horizontal turbulence and considered a general equation of state. They obtained 



U 2 (r) 



pgCpTiV^ - V + (<p/8)V,. 



(12) 



where M* = M (l — 2 ^ 0m ) * s the reduced mass and the other symbols have 
the same meaning as in Zahnl <1992ft and IMaeder fc Zahnl jl998). The driving 
term in the square brackets in the second member is Eq. It behaves mainly 
like E n ~ 



1 



2vrGp 



( gm ) • ^he term p means the average on the consid- 



ered equipotential. The term with the minus sign in the square bracket is the 
Gratton-Opik term, which becomes important in the outer layers when the local 
density is small. This term produces negative values of U 2 (r) (noted U(r) from 
now), meaning that the circulation is going down along the polar axis and up in 
the equatorial plane. This makes an outward transport of angular momentum, 
while a positive U(r) gives an inward transport. At lower Z, the Gratton-Opik 
term is neglig ible, which contr i butes to make larger J7-gradients in lower Z stars. 

Recently iMathis Sz Zahnl ( 2004f ) rederived the system of partial differential 
equations, which govern the transport of angular momentum, heat and chem- 
ical elements. They expand the departure from spherical symmetry to higher 
order and include explicitly the differential rotation in latitude, to first order. 
Boundary conditions for the surface and at the frontiers between radiative and 
convective zones are also explicitly given in this paper. 

In the ab ove equation D v is th e shear diffusion coefficient whose expression 
is taken as in iTalon fc Zaiml (|l997f l. The usual estimate of D h = i-r \2V(r) - 

aU(r)\ was given hv lZahnl (|1 992ft . Re cent studies suggest that this coeffi cient is 
at least an order of magnitude larger (|Maederll2003l : IMathis et afll2004ft . 



3. Equations for the transport of the chemical species 

The derivation of the equation for the transport of chemical species proceeds 
much like that of angular momentum: the mass fraction X{ of a given element 
i obeys an advection-diffusion equation which can be written 

^ + div(^U) = 
pD v r 2 ^ ) + ^-tt^ I sm6pD h ^- ) . (13) 



r 2 dr V V Or J r 2 sin 9 89 \ 89 
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Splitting the mass fractions in their horiz ontal average (Xj) and the ir fluctuation 
on the isobar, assuming that D h > D v , IChabover fc Zahnl (|1992T ) showed that 
Eq. EH can be written 

with the effective diffusivity 



D, 



<-ir 



rU (r) 



30D h 

D e Q expresses the effects of the meridional circulation and of the strong hori- 
zontal turbulence. This expression of D c q tells us that the vertical advection of 
chemical elements is inhibited by the strong horizontal turbulence characterized 
by Dh. For simplicity in Eg. 1141 the terms due to nuclear reactions are omitted, 
as well as the effects of element separation through radiation and gravitational 
settling. Let us finally add that various numeric al procedures f o r reso lving the 
diffusion equation are described and compared in lMevnet et al. 



4. Effects of rotation on the stellar winds 

The taking into account of the von Zeipel theorem in the frame of the radia- 
tive driyen w i nd theory has many inte resting consequences (jOwocki et al.lll996l : 
lMaederlll999l : iMaeder fc Mevnedl200oT ): 

1. The maximum luminosity of a rotating star is reduced by rotation. 

2. The expression for the critical velocity in a rotating star is different if the 
star is near or far from the Eddington limit. Far from the Eddington limit, 

■fcritjl = (firf) 5 ) where i? p b is the polar radius at break-up. When the 
Edd ington factor is bigg er than 0.639, the expression for t> C rit is different 
(|Maeder fc Mevnetl EoOO) . It is equal to 0.85, 0.69, 0.48, 0.35, 0.22, times 
v critt i for T = 0.70, 0.80, 0.90, 0.95, 0.98 and 1.00 respectively. 



Rotat ion induces anisotropi cs of the stellar winds ( Maeder fc Desiacauesl 
120011 : iDwarkadas fc Owoddll2002r t. 

Rotation increases the mass loss rates. We may estimate the mass loss 
rates of a rotating star compared to that of a non-rotating star at the 
same location in the HR diagram. The result is (|Maeder fc Mevnetl EoOO) 

£<2U (i-o*- , , (15 , 

M(0) r 1 i { y , 2 r i 

where T is the electron scattering opacity for a n on-rotating star wit h the 
same mass and luminosity, a is a force multiplier ( Lamers et al.l ll995). For 



B-stars far from the Eddington Limit, 4vr-r — 1.5 For stars close to T = 1 

to ' M(0) 

the increase of the mass loss rates may reach orders of magnitude. 
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It is interesting to note that three kinds of break-up limit can be defined de- 
pending on which mechanism, radiative acceleration, centrifugal acceleration or 
both contribute to counterbalance the gravity: 1— The T-Limit occurs when 
radiation effects largely dominate; 2.- The Q-Limit, when rotation effects are 
essentially determining break-up; 3.- The fir-Limit, when both rotation and 
radiation are important for the critical velocity. In the cases of Luminous blue 
variables the T- and Or-Limits likely occurs, while in the case of Be stars, the 
f2-Limit is the relevant limit. 



5. Some successes of rotating models 

The inclusion of rotation in massive star models improves in many respects 
the theoretical predictions. In particular, rotating models ca n reproduce the 
chem i cal enrichments observed at the surface of OBA stars ( Heeer fc Langeil 



2000|; iMevnet fc Maeder! 120001 V the number ratio of blue to red supergiants 
in the Small Magellanic Cloud (|Maeder Sz Mevnetll200lh . the variation with the 



metallicity of the number ratio of Wolf-Rayet to O type star s (IMevnet Sz Maedei 



2003) as well as of the type Ibc to type II supernovae ()Prantzos k, Boissiei 
200 31: IMevnet fc Maeder! 12004ft . Rotating models predict new chemical yields 
dHirschi et al.l 12004ft and naturally lead to the production of primary nitrogen 
at low metallicity (IMevnet Maede r! 12002ft . t hev also bring new views on the 
evolution of pop III stars ( Me vnet et al.ll2004aft . Comparisons with the observed 
rotational velocities of stars become possible and interesting consequences for 
the rotation rate of young pulsars and the progenitors of G amma Ray Burst can 
be obtained (|Wooslev fc Hegeril2004l : IMevnet et al.ll2004bft . 



6. Conditions for reaching break up during the MS phase 



Recently iTownsend et al. | (|2004h have given strong arguments indicating that 



the vsv&i measured for Be stars might be systematically underestimated. An 
important implication of their work is that these stars may be rotating much 
closer the critical velocity than generally assumed. As the authors note, rotation 
close to critical effectively reduces the effective equatorial gravity and could make 
material to easily leak into a disc. Typically, when v/v cr \t ~ 0.95, the velocity 
necessary to launch material into orbit is of the same order as the speed of sound 
in the outer layers. Therefor e mechanisms , as nonradial pulsation, begin to be 
effective for orbital ejection ( Owockill2004h . In view of these considerations, it 



appears interesting to investigate with our rotating stellar models, which are the 
initial conditions required for a star to reach the critical velocity during the MS 
phase. This is the object of this section. 

Figs. |21 show the evolution of the rotational velocities and of the fraction 
YY~ of the angular velocity to the critical angular velocity at the surface of star 
models of different initial masses between 9 and 120 Mq with account taken 
of anisotropic mass loss during the MS phase. The evolution of the rotation 
velocities at the stellar surface depends mainly on 2 factors, the internal coupling 
and the mass loss. 

1.- The coupling mechanisms transport angular momentum in the stellar 
interiors. The extreme case of strong coupling is the classical case of solid body 



Rotation and Mixing in Massive Stars 



9 




Age [10 6 years] Age [10 e years] 



Figure 2. Left /Evolution of the rotational velocities for star models of dif- 
ferent initial masses between 120 and 9 Mq with account taken of anisotropic 
mass loss during the MS phase. The rapid decrease near the end of the MS 
evolution is due to the bi-stability limit in the mass loss rates. The location 
of the end of the MS phase corresponds to the vertical decrease of the velocity. 
Right : Evolution of the fraction of the angular velocity to the critical an- 
gular velocity at the surface of star models of different initial masses between 
120 and 9 M© with account taken of anisotropic mass loss during the MS 
phase. 



rotation. In this case when mass loss is small, the star reaches the critical velocity 
during th e MS phase more or less q uickly de pending on th e initial rotation as 
shown by ISackman Anandl (|l970f ) (see also lLangerl I1997T ) . In the case of no 
coupling, i.e. of local conservation of the angular momentum, rotation becomes 
more and more subcritical. In the present models, with initial velocities at the 
beginning of the MS phase of ~ 260 km s -1 the situation is intermediate, with a 
moderate coupling due m ainly to meridional circulation, which is more efficient 
( Mevnet Maederl l2000l ) than shear transport, as far as transport of angular 
momentum is concerned. 

2.- For a given degree of coupling, the mass loss rates play a most critical 
role in the evolution of the surface rotati on. As shown by the co mparison of 
the models at Z = 0.02 and Z = 0.004 dMaeder fc MevnetJ Effll . for masses 
greater than ~40 Mq the models with solar composition have velocities that 
decrease rather rapidly, while at Z = 0.004 the velocities go up. Thus, for the 
most massive stars with moderate or strong coupling, the mass loss rates are the 
main factor influencing the evolution of rotation. The effect of mass loss is also 
well illustrated in Fig. El which shows the evolution of £l/Q CT i t at the surface of 
a 15 Mq solar metallicity model. One sees that during the MS phase, the ratio 
increases as a f unction of ti me, until logT e fj- becomes inferior to 4.40. According 
to Vink et al. (2000l: l200lh . below this value there is a sudden increase of the 
mass loss rates, which makes the rotation velocities rapidly decrease. 

It is interesting to note that below a mass of about 12 Mq, the mass loss 
rates are smaller and the internal coupling plays the main role in the evolution 
of the rotational velocities. This provides an interesting possibility of tests on 
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Figure 3. Three dimensional evolutionary track of a 15 M model at solar 
metallicity, with a velocity at the beginning of the MS phase of 260 km s _1 . 
The evolution in the theoretical HR diagram is plotted on the horizontal 
plane. The vertical axis shows the evolution of the surface angular velocity 
devided by the critical or break-up angular velocity (R c = Q/f2 cr it). 



the internal coupling by studying the differences in rotational velocities for stars 
at different distances of the ZAMS. In particular, such a study could allow us to 
test the role of magnetic coupling in radi ative envelop es, which is now a major 
open question in stellar rotation studies (|Sprurdl2002f ) . 

From Fig. |5J we see that the stars which have a velocity around 260 km 
s _1 at the beginning of the MS phase, do not reach the break-up limit. Taking 
into account of the wind anisotropies induced by rotation, a 9 M stellar model 
at solar metallicity with an equatorial velocity at the beginning of the core H- 
burning phase equal to 340 km s _1 reaches the break-up limit during the MS 
phase. This corresponds to values of O / O cr j t and of v/v CT n equal to 0.81 and 0.61 
respectively. The model with with an initial equatorial velocity of 290 km s _1 
(fi/f^crit = 0.71, w/fcrit = 0.52), approaches but does not reach the break-up 
limit during the MS phase. At the end of the MS phase, the 290 and 340 km s _1 
models present N/C ratios at the surface equal to 2.8 and 3.2 times the initial 
value. 

In Fig.^J evolutions oiv/v CT n for different initial metallicities and masses are 
shown. Interestingly, one notes that the 5 Mq stellar model at Z = 0.040 reaches 
the break-up limit well before the model at Z=0.020. This comes essentially from 
the dependence on density of the meridional velocity as given by the Gratton- 
Opik term. At higher metallicity, the radii of stars are bigger and thus the 
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Figure 4. Left : Evolution of the equatorial velocity divided by the crit- 
ical velocity at the surface of 5 M stellar models at different metallicities. 
RightT Same as in the left panel, for different initial mass stellar models at 
the metallicity Z = 0.001. The initial values of v q/v CI it are between 0.60 and 
0.64. 

density in the outer layers are smaller. At a given metallicity (see right panel 
of Fig. |1J), the higher initial mass star reach the break-up limit before the lower 
initial mass models. Again this can be explained through the density dependance 
of the meridional circulation velocity. 

Thus these numerical experiments show the following: 

1) For stars suffering little mass loss, the main effect which can bring these 
stars near the break-up limit is the outwards transport of angular momentum 
by the meridional circulation. Since the velocity of the meridional currents in 
the outer layers scales with the inverse of the density, the process becomes more 
efficient for stars of higher initial mass and/or higher initial metallicity. 

2) When the metallicity increases however, mass loss becomes more and 
more important and can prevent stars to reach the break-up limit. 

How the fraction of stars which reach the break-up limit varies as a func- 
tion of the initial metallicity requires a good knowledge of the distribution of the 
initial velocities as a function of the mass and of the metallicity. L arge and ex- 
tended surveys are being now performed with this purpose (see e.g. iNorth et alJ 
120041 : iRover et al.l 120041 ) . Also the study of the time evolution of the surface 
velocities on the Main Sequence for stars undergoing little mass loss would pro- 
vide very important hints on the efficiency of the internal tranport mechanisms 
of the angular momentum. 

Acknowledgments. We express our thanks to Raphael Hirschi for his help 
in producing Figure 01 
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